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This colloquium analyzes the interaction of light with two-dimensional periodic arrays of particles 
and holes. The enhanced optical transmission observed in the latter and the presence of surface 
modes in patterned metal surfaces are thoroughly discussed. A review of the most significant 
discoveries in this area is presented first. A simple tutorial model is then formulated to capture the 
essential physics involved in these phenomena, while allowing analytical derivations that provide 
deeper insight. Comparison with more elaborated calculations is offered as well. Finally, hole 
arrays in plasmon-supporting metals are compared to perforated perfect conductors, thus assessing 
the role of plasmons in these types of structures through analytical considerations. 
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I. INTRODUCTION 



The scattering of waves in periodic media plays a cen- 
tral role in areas of physics as diverse as low-energy elec- 
tron diffraction ( jPendry 1974) or atomic-beam scatter- 
ing from crystal surfaces (Farias and Rieder 19981. Va- 



lence electrons in solids, sound in certain ordered con- 



structions ( Martmez-Sala et al. 1995), or light in pho- 
tonic crystals ( Joannopoulos et al. 1997 Lopez 2003) 



limit their propagation in frequency regions known as 
band gaps (Ashcroft and Mermin 19761. Among these 



examples, the scattering of electromagnetic waves is par- 
ticularly important because it allows obtaining structural 
and spectroscopic information over a fantastically wide 
range of lengths, going from atomic dimensions in x-ray 



scattering (iHenke et al. 1993) to macroscopic distances 



in radio and microwaves. Actually, Maxwell's equations 
are written in first-order derivatives with respect to spa- 
tial coordinates, so that light scattering in the absence 
of nonlinear effects is solely controlled by the shape and 
permittivity of diffracting objects with distances mea- 
sured in units of the wavelength, and therefore, the same 
phenomena are encountered over entirely different length 
scales. 

We can classify the performance of periodic structures 
in three distinct categories according to the ratio of the 
period a to the wavelength A. For A ^> a, an effective ho- 
mogeneous medium description is possible. This is in fact 
what happens in most naturally-occurring substances 
when a has atomic dimensions. But also in certain arti- 
ficially textured materials (metamaterials), which allow 
achieving exotic behavior like magnetic response at visi- 
ble frequencies ( Grigorenko et al] |2005[) and media with 



negative refraction index ( Smith et al. 20041, without 



neglecting the exciting possibility of using nanoparticles 
as building blocks to tailor on-demand optical proper- 

The opposite limit (A <C a) is 
although 



ties ( Liz-Marzan 



_ 2006]. 

generally well accounted for by classical rays, 
keeping track of phases proves to be crucial near points 
of light accumulation, like in the self-imaging of grat- 



et al. 


2007 


Talbotj 


1836) 



undergo diffraction that under certain conditions can 



Nevertheless, it is the inter- 
mediate regime, when A is comparable to a, in which 
diffraction shows up in full display. We find examples 
of this in both three-dimensional (3D) photonic crys- 
tals, which offer a promising route to fully controlling 
light propagation over distances comparable to the wave- 
length (Joannopoulos et al. 1997 L6pez[[2003[ ), and two- 
dimensional (2D) crystals, in which an impressive degree 
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of optical confinement has been accomplished (Akahane derstand that it originates mainly in the interaction of 



et al 2003 1 



In this colloquium, we shall focus on light scattering 
by planar structures of particles or holes, which have 
become a current subject of intense research driven to 
some extent by advances in nano-patterning techniques. 
Our main purpose is to explain the phenomena observed 
within this context in a tutorial but nevertheless com- 
prehensive fashion. We shall first review experimental 
and theoretical developments in Sec. [TTJ Then, we shall 
formulate in Sec. |III| a simple powerful model that deals 
with the response of particle and hole arrays on a com- 
mon footing, leading to analytical expressions that cap- 
ture the main physical aspects of these systems. Finally, 
metals with plasmons will be discussed, and the main 
differences with respect to plasmon-free perfect conduc- 
tors elucidated, in Sec. 



IV 



We shall use Gaussian units, 
unless otherwise stated. 

The beginning of the last century witnessed impor- 
tant developments in diffraction of light in gratings af- 
ter Wood's observation of anomalous reflection bands 



(Wood 


1902, 


1912 


1935) 


tat ion 


Fano] 


1936 


1941 



1941| Lord Rayleigh 1907). Two 



types of anomalies were identified, one of them occurring 
when a diffracted beam becomes grazing to the plane 



of the grating, the Rayleigh condition (Lord Rayleigh 



1907), giving rise to a sharp bright band, and the other 
one showing up to the red of the former as an extended 
feature containing two neighboring dark and bright bands 



(Fano 1936 1941) 



The century concluded with another significant dis- 
covery (Ebbesen et al. 1998): periodic arrays of sub- 



wavelength holes drilled in thin metallic films can trans- 
mit much more light per hole at certain frequencies than 
what was previously expected for single openings, based 
upon Bethe's prediction of a severe cutoff in transmis- 
sion as (b/X) 4 for large A compared to the hole radius b 
(Bethe 1944). Previous knowledge gathered by electrical 



engineers in the microwave domain ([Chen. 1971 McPhe 



dran et al. 1980 Ulrich 1967) had already exploited the 



use of periodically-drilled surfaces as frequency-selective 
filters and discussed the occurrence of 100% transmis- 
sion at wavelengths slightly above the period. How- 
ever, the hole sizes that were considered in that con- 
text lied in the region of sizeable transmission for single 
holes. The more recently discovered extraordinary trans- 
mission phenomenon was however observed for narrower 
holes (relative to the wavelength), the transmission of 
which exceeded orders of magnitude what was expected 
from the sum of their individual transmissions ( |Ebbe-| 
sen et al. 19981. For square arrays under normal inci- 



dence, a transmission minimum occurred at a wavelength 
close to the period a, coinciding with the Rayleigh con- 
dition (Lord Rayleigh 19071, and a transmission maxi- 



mum showed up at longer wavelength, thus revealing its 
connection to Wood's anomalies 
|Sarrazin et al 



2003). However, t 



Ghaemi et al. 1998 



le explanation of the 



the apertures with surface plasmons ( 


Barnes et al 


2004 


Ghaemi et al, 1998 Martin-Moreno et al. 2001 


Popov 


et al. 


2000 


Salomon et al. 


2001 


Wannemacher 


2001), 



whereas other authors make emphasis in dynamical light 



diffraction ( Lezec and Thio 2004 Sarrazin et al. 2003 



Treacy |1999 2002). While the latter works well to un- 



derstand the observed extraordinary optical transmission 



in drilled plasmon-free perfect conductors ( 


Cao and Na- 


hatal |2004 Gomez-Rivas et aZ.[|2003||Mittra et aZ.[|1988| 


Miyamaru and Hangyo 


2004), supporters of the surface- 



plasmon interpretation argue that the enhanced trans- 
mission relies in this case on plasmon-like lattice-surface- 
bound modes sustained by p atterned perfect-conductor 
surfaces (|Pendry et al. 2004 1 . Actually, evidence of such 



modes had been observed before in periodically perfo- 
rated metallic screens for wavelengths several times larger 
than the period (Ulrich and Tacke 1972). We shall see 



below how these are in fact complementary views of the 
same phenomenon and how diffraction in particle arrays 
contains already the essential features that can be trans- 
lated to understand the phenomenology of hole arrays. 
But let us first summarize experimental and theoretical 
findings in this area. 



II. OVERVIEW OF EXISTING RESULTS 

A huge amount of literature has been accumulated on 
transmission through periodic structures, and it is an in- 
teresting exercise to reexamine it in connection to recent 
developments. 



A. Single holes 

Bethe's predicted cutoff in the transmission of a sin- 
gle hole in a perfect-conductor thin screen as (b/X) 4 is 
the leading-order term of the expansion of the transmis- 
sion cross section in powers of b/\ (Bel he 19441. Sub- 



sequent higher-order analytical corrections (Bouwkamp 



1954 Chang et al. 20061, and eventua l ly rigorous nu- 
meric al calculations ( [Garcia de Abajo 2002 Roberts 
19871, demonstrated that the cross section lies below 



the hole area up to a radius b m 0.2A. These results 
have found experimental corroboration down to the NIR 



regime (Obermiiller and Karrai 19951, with new local 



ized plasmon resonances showing up at shorter wave- 



lengths (Degiron et al. 2004 Rindzevicius et al. 2007) 



Two different mechanisms have been however sug- 
gested to achieve enhanced transmission in a single hole: 
filling it with a material of high permittivity (Garcia 
de Abajo[ |2002"[ |Garria-Vidal et al.\ [2005] |Webb and 



Li 20061, thus creating a partially-bound cavity mode 



that couples resonantly to incident light (see Sec. III.El; 
and decora ting the aperture with periodic corrugations 
(Lezec et al. 2002[) in much the same way as highly- 



effect is still a subject of debate, as some authors un- directional antennas are capable of focusing electromag- 
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netic radiation on a central dipole element by means 



of concentric, periodically-spaced metallic rings (James 



1977) 



B. Optical transmission through hole arrays 

The intensity of light passing through holes is boosted 
at certain wavelengths when we arrange them period- 
ically. Pioneering calculations and microwave experi- 
ments showed already zero reflection in thin films per- 
forated by periodic arrays of small apertures of radius 
b 



0.36A riChen 1971) 



Further seminal experiments 
focused on the relation between hole arrays in thin metal 
screens and their complementary screens (Ulrich 19671, 



putting Babinet's principle to a test in the far-infrared 
region. This was followed by numerous applied studies 
of hole arrays (regarded as frequency-selective surfaces) 
in the engineering community, including filters for solar 
energy collection and elements to enhance antennae per- 
formance flCwik et al] |1987| |Maystre[ |1980| |McPhedran 



et al. 1980 Mittra et al 



|1988| 



The work of Ebbesen et al. ( 1998 1 demonstrated in the 



optical domain extraordinary light transmission, which 
for the first time occurred for openings of radius below 
the cutoff of the first propagating mode in a circular 
waveguide, b < 0.29A. Since then, this phenomenon has 
been consistently observed for a varied list of metallic 



materials (Przybilla et al. 2006a I, over a wide range of 



wavelengths e.g., for microwaves ( Cao and Nahata 2004 



Gomez- Rivas et al. 20031, to which metals respond as 



nearly perfect conductors, in the infrared (Selcuk et al. 



2006), and in the VUV, using a good conductor in this 



regime like Al (Ekinci et al. 2007 1], and with different 
types of array symmetries, including the recent demon- 
stration of the effect in 2D quasi-crystal arrangements 



(Matsui et al. 2007 


Przybilla et al. 


2006b 


Schwanecke 


et al. 2006 


Sun et al. 20061. 



Two examples of enhanced transmission, taken from 
|Krishnan eTal] [2001] and |Martm-Moreno et al.\ pOOT] 
are shown in Fig. [1] The transmission is several times 
larger in the infrared peak than the prediction of Bcthc 
for non-interacting holes in a thin screen, and four or- 
ders of magnitude larger than what is expected for non- 
interaction apertures in a perfect-conductor film of the 
same thickness (dashed curves). 

Light transmission through hole arrays has been ex- 



amined theoretically for over four decades (Chen 1971 



Dawes et al\ [19891 [Eggimann and Collin| |1"962| |McPhe-| 



dran et al. 1980 1, although a detailed account of extraor- 



0.08 




Wavelength / Period Vf 

FIG. 1 (Color in online edition) Extraordinary optical trans- 
mission in hole arrays. The measured transmittance (solid 
curves) is shown for apertures drilled in gold (a) and silver (b) 
films, taken from |Kristman et al.\ |2001| and [Martin- Moreno 
\et al.\ |2001| respectively. The silver film is self-standing in 
air, while the gold is deposited in quartz and immersed in an 
index-matching liquid. The lattice constant is a — 600 nm in 
Ag and a = 750 nm in Au. The transmittance of the perfo- 
rated gold goes well above that predicted for non-interacting 
apertures in a perfect-conductor film or by Bethe's formula 
for a thin screen (dashed curves). Rayleigh's condition for 
the (1,0) and (1,1) beams becoming grazing are indicated by 
vertical dashed lines. Analytical results are shown as arrows 
in (a) and as a dashed curve in (b) (see Sec. IV. D I. The trans- 
mittance is presented vs the wavelength in the dielectric en- 
vironment of the metal, normalized to the lattice constant. 



the role of hole shape has been shown to yield nontrivial 



effects ( 


Elliott et al. 


2004 


Gordon et al] 


2004 Klein Ko- 


erkamp et al. 2004 


<rasavin et al. , 2005 van der Molen 


et al. 20051, such as larger enhancement and red shift 



of the transmission peaks with respect to the Rayleigh 
condition for light polarized along the short axis of elon- 
gated apertures. Finite arrays have been found to exhibit 
interesting shifts in the transmission maxima as well, de- 



dinary optical transmission in real metals had to wait un- pend i ng on the number of ap ertures ( |Bravo-Abad et al 



til the new century began ( Martin-Moreno et al. 2001 



Popov et al.\ [20 00 Salomo n et al.\ 2001| Sarrazi n et al 



2004 Lezec and Thio 2004). More exotic shapes like 



annular holes have been also simulated (Baida and Van 



2003 Wanncmachcr 2001) and the advance in compu- 



tation power allowed predictive capacity (|Chang et~al\ sured ( |Fan et aZ.[|2005p , with the additional appeal that 



2005 Klein Koerkamp et al. 2004 1 



Labeke 2002 Roberts and McPhedran 



1988) and mea- 



annular waveguides support always one guided mode at 



The influence of various geometrical and environmen- 
tal factors has been extensively studied. In particular, 



least (Jackson 19991 



The transmission is exponentially attenuated with hole 



depth because it is mediated by evanescent modes of 
the apertures regarded as narrow subwavelength waveg- 
uides. However, strong signatures of interaction between 



both metal interfaces have been reported ( Degiron et al. 
|2002[ ), as well as high sensitivity to dielectric environ 
ment, so that maximum transmission is achieved when 
the permittivity is the same on the two sides of the film 



(Krishnan et al. 2001) 



Extraordinary optical transmission has expanded to 



embrace a wide range of phenomena ( Genet and Ebbesen 



2007), like the interaction of hole arrays with molecules 



for potential applications in biosensing (Dintinger et al. 



2006a l and all - optica l switching (Dinting er et al.\ 2006b 
Janke et al. 2005 Smolyaninov et al. |2002 l, anc 
the demonstration of the quantum nature of plasmons 
through photon entanglement preservation after travers- 



ing a hole array ( Altewischer et al. 2002 1 



C. Particles 



The field of light scattering by small particles has a 



rich research tradition (Bohren and Huffman 1983 



de Hulst |1981h that is being continued by hot topics 



such as for example novel near-field effects in the coupling 
of metallic nanoparticle arrays (Krenn et al. 19991 and 
strong inter-particle interactions in dimers (|Atay et aT\ 



|2004||Nordlander et q?.||2004||Romero et al 



200G). Here, 



we shall single out just two recent exciting developments 
in line with the rest of our discussion. The first one refers 
to coupled metallic nanoparticle arrays. These particles 
can sustain localized plasmon excitations that hop across 
neighbors. It has been suggested (Quinten et al. [19981, 



and later confirmed by experiment (Ma ier et aL[ 2001 
|2003[ ), that this phenomenon can be utilized to transmit 
light energy along chains of subwavelength particles, thus 
providing some basic constituents for future plasmonic 
devices. 

In a different development, the scattering spectra from 
ID and 2D arrays of metallic nanoparticles were pre- 
dicted to exhibit very narrow plasmon lineshapes pro- 
duced by d y namic al scattering (Zou et al. 2004| [Zou 



and Schatz 20041. Experiments performed on litho- 



graphically patterned particle arrays confirmed this effect 
and achieved reaso nable control over spectral lineshapes 
We shall discuss this further in Sec. 



2005) 



(Hicks et al 
IIII.A.2I 



III. TUTORIAL APPROACH 




E ex, (R) 



y * 



* * K * 



FIG. 2 (Color in online edition) Two-dimensional array of 
small identical particles illuminated by a light plane wave. 
k|| is the momentum component parallel to the array. The 
particle at position R n displays a dipole p n . 



and making several challenging predictions. 



A. Basic relations 

We shall start with some basic analytical relations for 
the scattering of an external light plane wave on a peri- 
odic array of identical particles that are small compared 
to both the wavelength and their separation (see Fig. 
[2|. Within linear, non-magnetic response, the particle 
at position R„ can be assumed to respond with an in- 
duced dipole p„ = oeE(R„), determined by its electric 
polarizability tensor a~E and the self-consistent field act- 
ing on it, E(R„). This dipole induces an electric field at 
point r that can be written Q°(r — R„)p„ in terms of the 
dipolc-dipolc interaction tensor, 



g°(r) = (fc 2 + VV)- 



Akr 



(i) 



where k is the light momentum in free space. 1 Now, the 
self-consistent dipole of our particle is found to be 



Pn = "E 



E cxt (R„)+ £°(Rn-R„'W 



, (2) 



where E cxt (R„) = E cxt exp(ik|| -R„) is the external elec- 
tric field, which depends upon the site position R n just 
through a phase factor involving components of the in- 
coming wave momentum parallel to the array, ky , as illus- 
trated in Fig. [2] and the second term inside the square 
brackets represents the field induced by the rest of the 
particles. Bloch's theorem guarantees that the solution 
of Eq. (|2j) must have the form p„ = pexp(iky • R„). 
Direct insertion of this expression into Eq. ([2| leads to 



A tutorial model will be presented next that becomes 
exact in the limit of narrow holes or small particles in 
perfect-conductor films. This model will describe the ba- 
sic physics involved both in extraordinary light transmis- 
sion and in lattice surface modes of structured metals, 
but it leads to simple analytical expressions that permit 
understanding these phenomena in a fundamental way 



l/a E -G(k,|)" 



-E c 



(3) 



1 More explicitly, Q°(r)p = [exp(ifcr)/r 3 ] { [(fcr) 2 + ikr - l] 
[(fcr) 2 + 3ifcr - 3] (r-p)r/r 2 }. 
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and 



G(k|, 



^S (R„K 



ik M R„ 



(4) 



where we have chosen Ro = 0. Notice that the denomi- 
nator of Eq. ([3]) separates the properties of the particles 
(cue) from those of the lattice [the structure-factor-type 
of sum G(k||)], in the spirit of the KKR method in solid 
state physics ( Ashcroft and Mermin[ 19761. The lattice 
sum in Eq. Q can be converted into rapidly converging 



sums using Ewald's method (Glasser and Zucker 1980), 



and we have used in particular the procedure elaborated 
by |Kambe|fll968[ ). 

Incidentally, Eqs. Q-Q can be also applied to 3D 
particle arrays with kii replaced by a 3D crystal momen- 
tum. This type of approach has been shown to lead to 



robus t band gaps in atomic lattices ( van Coevorden et al. 
1996). Furtherm ore Eq. |2| together wi th the Clausius 



Mossotti formula ( jAshcroft and Mermin 1976 1 constitute 
the basis of the discrete-dipole approximation (DDA) 
method for solving Maxwell's equations in arbitrary ge- 



ometries (Draine and Flatau 



1994 Purcell and Penny- 



packer 19731. It should also be noted that the present 
approach can be extended to larger particles arranged in 



ordered (Stefanou et al. 1998 2000 1 or disordered arrays 



( Garcia de Abajo 1999 1 by including higher-order mul- 



tipoles, and that this is one of the methods that can be 
actually applied to d educe effective optical properties of 
composite materials ((Milton 20021. 



It is useful to represent the dipole-dipole interaction in 
2D momentum space in the plane of the array, which we 
shall take to coincide with z = 0. This is readily done by 
expressing the scalar interaction at the right end of Eq. 
@ as 



Akr 



2tt 



d 2 Q 



,i(Q-R+ft*|z|) 



where k z = \Jk? — Q 2 is the normal momentum and the 
notation r = (R, z), with R = (x,y), has been adopted. 
From here and Eq. Q one obtains expressions like 

2n / k 



GL(r) 



2\ i(Q-R+fe 2 |z|) 



(5) 



for the components of the interaction tensor, here spec- 
ified for the xx directions. This allows us to recast Eq. 
Q into a sum over 2D reciprocal lattice vectors g, using 
the relation 

]Texp(iQ.R n ) = ^]r<5(Q-g), ( 6 ) 

n g 

where A is the area of the lattice unit cell. For example, 
the G xx component under normal incidence {k\\ = 0) 
becomes 



G xx (0) = lim 



1 

2^ 



K z 



A 
t 

d 2 Q 



(fc 2 -Q2 )e ifc 2 M 



where fcf = y/k 2 — g 2 , and the integral represents the 
subtraction of the n = term in the sum of @. This 
expression is important to elucidate some properties of 
the lattice sums, as we shall show below. 



1. Reflection and absorption in particle arrays 

The scattered field is given by a Rayleigh expansion 
simila r to the one in Eq. Q ( jGarcia de Abajo et al. 



2006 1 , with each vector g labeling one reflected and one 
transmitted be am of parallel momentum kn + g (Lord 
Rayleigh 19071. In the far field in particular, the zero- 



order (g = 0) reflection and transmission coefficients un- 
der normal incidence reduce to 



2nik/A 



and 



1/oe-Gx«(0) 



t = 1 



(8) 



(9) 



where the first term in the right-hand side of Eq. ^ 
represents the unscattered beam, and the numerator of 
^ is the far-field amplitude produced by a lattice of unit 
dipoles. 

Interestingly, the absorbance of the array is given by 
1 — |1 + r\ 2 — \r\ 2 [see Eq. (pj], which when regarded as 
a function of the complex variable r, has a maximum of 
50% coinciding with r = —1/2 and t = 1/2. This con- 
dition is easily attainable near a lattice singularity (see 
Sec. III.B I, using for instance weakly dissipative spherical 



particles. Similar results have been predicted for narrow 
cylinder arrays (Laroche et al. 20061, in which 100% ab- 



sorption is possible in one of the polarization components 
for the right choice of parameters. 

A particularly simple situation is encountered when 
the wavelength is larger than the lattice spacing, so that 
all diffracted beams other than the zero-order beam are 
evanescent (|k|| + g| > k). Then, upon inspection of Eq. 
([7]), one finds the useful relation 



3{G ra (0)} = 2irk/A- 2fc 3 /3, k < g u 



(10) 



where g\ denotes the period of the reciprocal lattice 
(<7i = 2n/a for square arrays). Moreover, if the par- 
ticles are non-absorbing, the optical theorem constrains 
their polarizability by the condition 3{-1/oe} = 2A ,3 /3 
(van de Hulst 1981). Combining these expressions, one 



obtains 



-1 



iA 
2-rrfc 



5J{l/a E -G ra (0)} 



(11) 



for the reflection coefficient of non-dissipative particles 
under normal incidence below the diffraction threshold. 

The electrostatic approximation provides a reasonable 
description of the electric polarizability of small parti- 
cles, dig. However, this needs to be amended in order 



6 




1.00 



1.02 



1.04 



1.06 




1.0000 1.0005 1.0010 1.0015 
Wavelength / Period 



1.0020 



FIG. 3 (Color in online edition) Reflectance spectra of square 
arrays of perfectly-conducting thin circular disks. The wave- 
length A is normalized to the lattice constant a. The disks 
radius is b = a/5 in (a) and 6 = a/9 in (b). The light is im- 
pinging normal to the array and 100% reflection is observed in 
these two cases at the maximum. Solid curves: full nu mer ical 
results. Dashed curves: analytical model for |r| 2 [Eq. (111]. 



to comply with the mentioned optical-theorem constrain, 
for instance via the prescription a-E — — 2ifc 3 /3). 

Analytical expressions for q| ! exist for a variety of parti- 
cle shapes, including homogeneous spheres (a| s = 6 3 (e — 
l)/(e + 2), where b is the radius and e the permittivity) 
and ellipsoids (Jones 1945). 



We illustrate the applicability of Eq. (11) through 
an example consisting of square lattices of perfectly- 
conducting thin disks. Fig. [3] compares the analytical 



result of Eq. (11) (dashed curves) with the full solution 
of Maxwell's equations obtained by following a layer- 
KKR multiple-scattering formalism (Stefanou et al. 
fl998| , |2000| to simulate the array together with a modal 
expansion solution of the isolated disk similar to the one 



available for isolated holes (Garcia de Abajo et al. 
|Roberts| 



2005a 



1987). In the analytical solution we have used 



the polarizability of thin metallic disks as derived from 
an ellipsoid of vanishing height, a| s = 46 3 /37r, where b is 
the radius. The results of the analytical model describe 
qualitatively well the presence of zero- and full-reflection 
points in the spectra, irrespectively of the disk size, but 
we shall discuss this point further in Sec. III.C 



2. Narrowing lineshapes through dynamical scattering 

The above formalism can be used to explain the effect 
of narrowed plasmon lineshapes in the scattering spectra 



of ID and 2D particle arrays (Hicks et al. 2005 Zou 



et al. 2004 Zou and Schatz 20041. For simplicity, we 



shall discuss metallic: spherical particles described by the 
Drude dielectric function 



Uj(uO + IT}) 



(12) 



where uj p is the bulk plasma frequency and the plasmon 
damping rate is w rj/2 -C oj p . 

Using this expression to obtain the polarizability of 
a small sphere of radius 6 (see Sec. III. A.l I, we can 
recast Eq. (J3j> into a Lorentzian of width w r}/2 + 
(uj p b 3 /2v / 3)S{G} . The natural width of the isolated par- 
ticles is now supplemented by a term proportional to 
3{G} [see Eq. (JToJ)] , which can take negative values that 
compensate the r}/2 term to render arbitrarily narrow 
collective plasmon resonances for an appropriate choice 
of array parameters. 

Applying this to a 2D square array under normal inci- 
dence with A ~ a, we find that Eq. ( 10 1 yields complete 
cancelation of the width for b/a « .16(77 / ) 1 / 3 . Under 
such conditions, the narrowing of the width is just lim- 
ited by the physical requirement that \r\ 2 + \t\ 2 < 1 [see 
Eqs. (j§ and 



B. Lattice singularities 

The interaction among particles in the periodic arrays 
of Sec. |III.A| appears to be governed by the lattice sums 
G(k|| ) and is dominated by their singularities, which orig- 
inate in accumulation of in-phase scattered fields. Fol- 
lowing similar arguments to previous expositions of this 
idea (Fano 1941| Lord Raylcigh 19071, we just consider 



a ID periodic chain of particles illuminated by an in- 
cident plane wave with both propagation direction and 
electric field perpendicular to the array, so that the field 
induced by a given particle on a distant one scales with 
the inverse of their separation, and thus, the contribution 
of distant particles to the interaction lattice sum has the 
convergence properties of the series Y^=i e lkan /n, which 
diverges as th e wavelength approaches the period a as 
— ln(fca — 27t) ( Gradshteyn and Ryzhik 1980). The same 



is true for 2D arrays. These singularities in G(kii) arc 
signaled by the Rayleigh condition of a diffracted beam 
becoming grazing ( Lord Rayleigh[ 1907), as can be seen 
from Eq. ([7]), where divergent terms g rs k (i.e., terms 
with zero normal momentum k 9 z ) dominate the sum. 

A remarkable consequence of this analysis is that the 
array becomes invisible to the incoming light right at the 
lattice sum divergence \G XX (0) — > 00, so r — > 0, according 
to Eq. (|8j], showing 100% transmission even for absorb- 
ing particles. 
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FIG. 4 (Color in online edition) Lattice sum G zz (k||) [Eq. 
Q] for a square lattice of period a as a function of parallel 
momentum fcii and wavelength A. The direction of kn is along 
one of the axes of the lattice. 



Focusing for simplicity on a square array of period a, 
the normal-incidence la ttice sum ^ diverges as (Garcia 
de Abajo et a?.[|2005a| 



G xx (0) 



V x / a - 1 



118 



(13) 



for A > a, where a fitted constant has been subtracted in 
order to extend the validity of this expression well beyond 
the singularity. 

For oblique incidence with kn along one of the lattice 
unit vectors, proceeding as in the derivation of Eq. ([7]), 
one finds that G(kii ) is diagonal and its components di- 
verge as 



1 



+ 2im/a) 2 + {2irl/a) 2 - W ' 



(14) 



where n and I run over integral numbers (excluding I = 
in G xx ). This behavior is illustrated in Fig. |1J show- 
ing in full display the lattice singularities exhibited by 
^{G Z M}. 



C. Hole arrays 

1. Babinet's principle and hole arrays in thin screens 

The behavior of hole arrays in perfect-conductor 
screens can be directly connected to the properties of the 
disk arrays considered in Fi g. [3) Indeed, one can invoke 
the exact Babinet principle ( Born and Wolf[ 1999 Jack- 



son 



1999 ) , which connects the reflected fields of the disk 



array for a given incident polarization with the transmit- 
ted fields of its complementary hole array with orthogo- 
nal polarizati on, as illustrated in Fig. [5] (Garcia de Abajo 
et al.\ 2005a I. Therefore, the reflectance spectra shown 
in Fig. 3 arc identical with the transmittance spectra of 
the complementary perforated screens. 




1.00 1 .02 1.04 

Wavelength / Period 



1.06 



FIG. 5 (Color in online edition) Babinet's principle applied 
to disk and hole arrays. The transmittance (reflectance) of the 
disk array for light of a given polarization a (s or p) is identical 
to the reflectance (transmittance) of the complementary hole 
array for orthogonal polarization a' (p or s, respectively). 



Focusing again on square arrays and normal incidence, 
we observe two characteristic features in the transmit- 
tance spectra: (i) the transmission vanishes when the 
wavelength A equals the period a, and (ii) a 100% trans- 
mission maximum takes place at a wavelength slightly 
above o. The origin of these effects can be traced back to 



Wood's anomalies in gratings (Wood 1902 1912 19351 



and to their interpretation in terms of the following two 
mechanisms (Fano 1936 19411: (i) accumulation of in- 
phase scattering events when the wavelength equals the 
period (see explanation in Sec. |III.B I and (ii) coupling 
of the incident light to a surface resonance. These phe- 
nomena persist in hole arrays perforated in thicker films 
of non-ideal absorbing metals, for which the maximum 
transmission is reduced but still justifies the term ex- 



traordinary optical transmission (Ebbesen et al. 1998) 



The analytical simplicity of the transmission coefficient 
for our thin-screen hole array, given by the right hand side 
of Eq. (Ill, allows us to gain deeper insight into the origin 
of this phenomenon. The lattice sum G xx (0) was shown 
to diverge when A — a, as Fig. [^illustrates. This leads to 
vanishing transmission, which we can interpret in terms 
of accumulation of in-phase scattering (see discussion in 
Sec. III.B I. Furthermore, 100% transmission is achieved if 



the second term in the denominator of Eq. (Ill becomes 



zero, a condition that can be rigorously fulfilled for ar- 
bitrarily tiny apertures ( Garcia de Abajo et al.\ 2005a I: 
the smaller the holes, the larger TJoe, because the polar- 
izability is proportional to the cube of their radius, but 
no matter how large this fraction becomes, there is al- 
ways one wavelength at which the divergent lattice sum 
matches it. This statement is illustrated by geometrical 
construction in Fig. [6j in which the point of intersection 
of the horizontal dotted line and the solid curve [Fig.|6ja)] 
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respond to Fano profiles of the form (Fano 1961 1 



i di.sk ! I 



t holes 



-» T. 



holes 



FIG. 6 (Color in online edition) Geometrical construction of 
the condition of full transmission in a hole array, (a) Wave- 
length dependence of the real part of the lattice sum G xx [Eq. 
Q] for fey = 0. (b) Normal-incidence transmittance of a hole 
array complementary of the disk array of Fig.[3j[a) (b = a/5): 
exact calculation (solid curve), analytical model of Eq. (Ill 
(dashed curve), and Fano profile of Eq. ( |15[ ) (dotted curve). 
The transmission minimum at A = a results from the diver- 
gence of Gn, while the transmission maximum (see vertical 
dashed line) is derived from the condition that $1{G XX } equals 
the inverse of the hole polarizability, according to Eq. (111. 



signals the condition ^{l/a^ — G xx (0)} — 0. 2 The pos- 
sibility of 100% transmission in non-absorbing structures 
has been pointed out before (MaystreJ 1980 McPhedran 



et al. 19801), and the theory just presented goes further to 
show that this is possible for arbitrarily small holes. Nev- 
ertheless, the number of apertures needed to accomplish 
high transmission will increase as they become smaller, 
and at the same time the transmission resonance will 
be increasingly narrower and closer to X — a. There- 
fore, these transmission maxima involve long-range in- 
teraction among holes, dominated by dynamical diffrac- 
tion (i.e., multiple-scattering paths). In fact, if only 
single-scattering were considered, Eq. (J3j> would become 
p = oe (1 + aEG(k||)aE)E cxt , which wrongly predicts si- 
multaneous divergence of transmittance and reflectance 
at A = a. 

This collective response in planar periodic arrays can 
be regarded as a lattice surface resonance (Fano 19411, 
which becomes a true surface-bound state when evanes- 
cent incoming waves are considered, as we shall see in 



Sec. III.D However, the resonance is strongly coupled to 
propagating light for external plane-wave illumination, a 
situation described by Fano (Fano 1961 1 in his study of a 



discrete resonance state (our lattice surface-bound mode) 
coupled to a continuum (the transmitted light). This 
type of approach has been shown to work rather well in 
theory ( Chang et al. 2005| |Sarrazin et al. 2003[ ) and in 
comparison with measured transmission spectra (Genet 
et al. 20031. Our transmittance calculations should also 



T = C 



{q + ef 



(15) 



where e can be assimilated to the light frequency and q 
describes the strength of the coupling to the lattice sur- 
face resonance. Fig. [6lb) compares our exact calculation 
of the transmittance (solid curve) with a Fano profile cor- 
responding to parameters q — —3 and C = 0.1 (dotted 
curve) , in which we assume a linear relationship between 
e and the light frequency, with e = —0.33 (e = 3) for 
T = 1 (T = 0). The agreement is very reasonable, con- 
sidering that no dependence of the coupling parameter 
on wavelength is taken into account. This further sup- 
ports an interpretation of extraordinary transmission in 
terms of coupling to the lattice surface resonance set up 
by dynamical diffraction in the array. 

The geometrical construction of Fig. [6] provides a vi- 
sual explanation of transmission in arrays of elongated 
apertures: an elongated piece of planar metal (e.g., a 
rectangle) has larger electric polarizability along its long- 
axis direction, and this has direct consequences for the 
Babinct-rclatcd situation of an elongated hole with the 
electric field along the short axis; larger polarizability 
involves more red-shifted and broader transmission max- 
ima [this is so because the point of crossing in Fig. [6^ a) 
occurs wh ere G xx is less steep] , just as observed exper 



iment ally \ Gordon et al. 
20041. 



2004 Klein Koerkamp et al. 



Incidentally, Eqs. ([3| and Q constitute a good ap- 
proximation to describe the extraordinary transmission 



observed in 2D quasi-crystal hole arrays (Matsui et al. 
2007] |Przybilla et al.\ |2006b| |Schwanecke et al 



2006 



Sun et al. 2006), in which the lattice sum G exhibits 



pronounced, but finite maxima related to bright spots 
in the Fourier transform of the hole distribution. These 
spots define the reciprocal lattice for periodic arrays, but 
have quasi-crystal angular symmetry in quasi-crystals. 
In the spirit of Raylcigh's explanation of Wood's anoma- 
lies (Lord Raylcigh, 1907), the cumulative effect of long- 



distance interaction among apertures can be claimed to 
create these reciprocal-space hot spots, so that the ef- 
fect of neighboring holes can be overlooked and an ef- 
fective homogeneous p describes qualitatively well the 
extraordinary transmission effect in quasi-crystal arrays 
(Schwanecke et al. 20061, as well as the rich Talbot- like 



structure and subwavclcngth light localization observed 
at distances up to several wavelengths away from the ar- 
ray (Huang et al. 20071. 



2. Single holes in thick films 



2 We rely here on the condition 5R{1/«e} > 0, which is satisfied 
by the polarizability of planar, perfectly-conducting disks. Inter- 
estingly, lattice resonances will be absent in arrays of particles 
with negative polarizability, such as metallic nanoparticles under 
blue-detuned illumination relative to a nearby plasmon band. 



Our use of Babinet's principle in the previous section 
indicates that, similar to small particles, small holes in 
perfect conductors can be assimilated to equivalent in- 
duced dipoles, in line with Bethe's pioneering description 
of the field scattered by a single aperture in a thin screen 
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(Bethe 19441, which he regarded as arising from a mag- 
netic dipole parallel to the screen plus an electric dipole 
perpendicular to it. 

Narrow holes can still be represented by induced 
dipoles in thick screens, as illustrated in Fig. [7]ja). Par- 
allel electric dipoles and perpendicular magnetic dipoles 
are forbidden by the condition that the parallel electric 
field and the perpendicular magnetic field vanish at a 
perfect-conductor surface. This allows defining electric 
(E) and magnetic (M) polarizabilities both on the same 
side as the applied field (a„, with v =E,M) and on the op- 
posite side (o/ u ). Furthermore, energy flux conservation 
under arbitrary illumination leads to an exact optical- 
theorem type of relationship between these polarizabili- 



ties (Garcia de Abajo et al. 2006): by considering two 



plane waves incident on either side of the film and by im- 
posing that the incoming energy flux equals the outgoing 
one (because perfect conductors cannot absorb energy), 
we obtain the condition 



-2k 3 



(16) 



where we have defined 



±a' 



as hole response functions. The remaining real parts of 
are obtained numerically from the field scattered by 



a single hole (Garcia de Abajo 2002 Roberts 1987). 



These functions are represented in Fig. [7[b)-(c) within 
the electrostatic limit, clearly showin g |?a{ff^i~} | ~~ * 00 m 
the thin film limit, where a' v = — a v (Jackson 19991. 



3. Hole arrays in thick films 

Periodic arrays of sufficiently narrow and spaced holes 
can also be described by perpendicular electric dipoles p 
and p' and parallel magnetic dipoles m and mf, where 
primed (unprimcd) quantities are defined on the en- 
try (exit) side of the film, as determined by the in- 
coming light [see Fig. |7fa)]. We consider first a unit- 
electric-field p-polarized plane wave incident on a hole 
array with parallel momentum ky along x, so that the 
external field (incident plus reflected) in the absence 
of the apertures has parallel magnetic field H° xt = 2 
along the y direction and perpendicular electric field 
E° xt — —2kn/k along z. Then, one can generalize Eq. 
([3]) and write a set of mul tiple-scattering equations for 



the self-consistent dipoles (Collin and Eggimann 1961 
|Eggimann and Col lin 1962). Symmetry considerations 



demand that our magnetic and electric dipoles be ori- 
ented as m = my and p = pi, respectively. Following 
the notation of Sec. |III.A[ we can write 
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FIG. 7 (Color in online edition) Response of a small hole 
in a perfect-conductor thick film, (a) The field scattered 
by a subwavelength aperture in response to external electric 
(E ) and magnetic (_ff oxt ) fields is equivalent (at large dis- 
tance compared to the radius 6) to that of effective electric 
(p) and magnetic (m) dipoles, which allow defining polariz- 
abilities (ftE and cim, respectively) both on the same side as 
the external fields (a v ) and on the opposite side (a'„). Only 
the perpendicular component of the electric field and the par- 
allel component of the magnetic field can be nonzero at the 
surfaces of the perfect-conductor film, (b)-(c) Thickness de- 
pendence of the real part of the hole response functions 
for A > 6 [the imaginary part satisfies Eq. 1 16 (]. 



with a new lattice sum defined as 

H = -ikJ2^~ ikllXn d Xn 



This sum stands for the interaction between mixed elec- 
tric and magnetic dipoles. We can understand the above 
equations in a very intuitive way; for instance, the first 
one of them states that the electric dipole on the entry 
side (p) results from the response to the z-component of 
the self-consistent field on that side (E° xt + G zz p — Hm) 
via the polarizability cve plus the response to the self- 
consistent field on the opposite film side (G zz p' — Hm') 
via Qg. The solution to these equations can be readily 
written as 



p±p' = -2[(<?± 
m±m' = 2[( 5 ± - 



-G yy )k {] /k + H]/A±, (17) 
G zz ) + HkJk]/A±, (18) 
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with 



H 2 



A± = (<7e - G zz )(g M - G V y) 

The zero-order transmittance of the holey film is then 
obtained from the far field set up by the infinite 2D ar- 
ray of induced dipoles, T p = \(2irk 2 / Ak z )(m' — p'k\\/k)\ 2 , 

where k z = * k 2 — k 2 . 



Similar considerations for s-polarized light show that 
its transmittance reduces to T — \2-Kkm' / A\ 2 , with mag- 
netic dipoles parallel to k|| and no electric dipoles whatso- 
ever (E° xt — 0). More precisely, m±m' = (2k z /k)/(g^ l — 
G xx ), from which one obtains 



1 



2?rfc 2 



2-nk 



9u 
1 



— Gn 



9m 



(19) 



G X x} 1 



2ivk z 



1 



for the transmittance. The last identity in Eq. ( 19 I comes 
from Eqs. (10 1 and (16) for diffractionless arrays. 



Interestingly, Eq. (19) predicts 100% transmission if 



1 



A 

2nk, 



- G xx } M{g M - Gxx} = 0. (20) 



This is a second-order algebraic equation in %i{G x 
admits positive real solutions provided 

+ _ -I 
Ank \9m 9m\ 



> 1. 



,} that 



(21) 



Actually, ^i{G X x} can match those roots near the I ^ 
singularities of Eq. ( 14 1 , where it can be chosen arbitrar- 



i ly l arge within a narrow range of wavelengths [see Eq. 
( 13 >] . It should be noted that the difference g^j — 
falls off rapidly to zero when the film thickness h is made 
much larger than the hole radius b [see Fig.^b)]. How- 
ever, if we fix both the h/b ratio and the angle of in- 
cidence, the left hand side of (21 1 reduces to a positive 
real constant times XA/b 3 , leading to the conclusion that 
100% transmission is attainable at a wavelength close to 
the Rayleigh condition (e.g., A > a for normal incidence 
on a square lattice of spacing a) regardless how narrow 
the holes are as compared to the film thickness. Surpris- 
ingly, this requires that the ratio of the lattice constant 
to the hole radius be increased for deeper holes in order 
to compensate the fall in g~j^ — g~^ for larger h/b. 

The transmittance shows an interesting dependence on 
film thickness h ( |Martfn-Moreno et al. 2001[ ), as illus- 
trated in Fig. [8j The maximum of Fig. [6j is initially 
blue-shifted closer to A — a for small h, accompanied 
by a second narrower peak at even s hor ter wavelengths 
3 [these are the two solutions of Eq. ( 20 ) under the con- 
dition (21 1]. As h increases, inter-side interaction weak- 
ens and the two 100% maxima approach each other. At 
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FIG. 8 (Color in online edition) Thickness dependence of 
the normal-incidence transmittance spectra of square arrays 
of circu lar holes drilled in perfect-conductor films, according 
to Eq. ( 19 1. The hole radius b = 0.2a, the wavelength A, and 



the film thickness h are given relative to the period a (see text 
insets) . 



some point only one transmission maximum is observed 
when the left hand side of (21) is exactly 1. For even 
thicker films, the condition (21 ) cannot be met any longer 

The 



and the transmission maximum departs from 100%. 
Fano character of these lattice resonances is again visible 
through vanishing transmission at a wavelength immedi- 
ately below the maximum. 

Incidentally, perfect conductors are perfectly non-lossy, 
so that light dissipation must take place only at the open- 
ings if they are infiltrated with some dissipative material. 
For deep enough holes, the transmission is negligible and 
the absorbance becomes 1 — |r| 2 , which can reach 100% 
values under suitable resonant conditions, for instance in 
the IR by combining holes drilled in noble metals (behav- 
ing nearly as perfect conductors) infiltrated with phonon- 
polariton materials. In fact, a similar effe ct has been 



observed in th e visible using Au gratings (Hutley and 
|Maystre| |1976[) and in the infrared using SiC gratings 



OGreffet et al. 



2002J). 



D. Lattice surface modes in structured metals 

The flourishing area of plasmonics is demonstrating 
how confining electromagnetic fields to a surface can find 



many potential applications on the nanoscale (Ozbay 



2006[). Z enneck waves at radio frequencies (Barlow 1958 
1907), phonon-polaritons in the infrared (Gr 



Zenneck 



effet et al. 2002 Hillenbrand et al. 2002), and plas 



language of Fano arise from coupling to different light continua on 
either side of the film, but one of these resonances has vanishing 
3 In fact, there are two lattice resonances for h = 0, which in the width and is placed at A = a due to strong inter-side interaction. 
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mons in the visible are in fact different manifestations 
of the same phenomenon: confinement of electromag- 
netic fields to curved or planar surfaces. Even perfect- 
conductor screens, which are unable to trap light when 
they are flat, were experimentally shown by |Ulrich and| 



Tacke (19721 to host confined surface modes of p polar- 
ization when molded into films pierced by periodic arrays 
of holes spaced a distance much smaller than the wave- 
length [see Fig. ^b)]. 

In a recent independent development, |Pendry et al.\ 
( 2004 ) have studied surface modes in drilled semi- infinite 



metal, suggesting the possibility to extend plasmon-likc 
behavior to lower-frequency domains via the flattening 
of the mode dispersion relation driven by propagating 
modes of the holes, and stimulating new microwaves 

The analysis of 



observations 


Hibbins et al. 


Pendry et al. 


(20041 relied on 


based upon t : 


leir lowest-order 



2005) 



modes), which allowed extracting local permittivity and 
permeability functions in a metamaterial approach to ho- 



ley metals. However, |Garcfa de Abajo and Saenz (20051 
showed later that higher-order modes (and in particular 
TM modes) are important, giving rise to large quantita- 
tive modifications to the dispersion relation and revealing 
finer details in the holey metal response that go beyond a 
simple local metamaterial description (e.g., the angular 
dependence of the reflection coefficient does not follow 
the Fresnel equations with local optical constants). 

At variance with planar perfect conductors and their 
lack of surface modes, corrugated metallic surfaces can 
support bound states even in the long- wavelength limit. 
In an intuitive picture, surface confinement in a drilled 
semi-infinite perfect conductor can be related to the 
evanescent penetration of the electromagnetic field in- 
side the holes, in much the same way as surface plasmons 
enter a distance of the order of the skin depth inside a 



metal in th e visible and NIR regimes (Barnes and Sam- 
|bles| [2004 ) . Actually, these modes share with plasmons 
their character of p-polarized evanescent waves. 

Next, we elaborate a tutorial, analytical formulation of 
this phenomenon that becomes exact in the limit of small 
holes of size s <C a < A, arranged in a lattice of period a 
(Garcia de Abajo and Saenz 20051. Although we focus 



our analysis on periodic hole arrays drilled in a semi- 
infinite perfect-conductor, it must be emphasized that 
periodicity is not really needed and that similar modes 
should exist for patterns other than holes (e.g., small 
protuberances or particles deposited on an otherwise flat 
surface). 

Using the formalism of Sec. |III.C[ we find that Eqs. 



( 17 1 and ( 18 1 offer a simple description of lattice surface- 
bound modes in metallic films. For infinitely-deep square 
holes as sketched in an inset of Fig. [9]ja), the surface 
modes must correspond to non-vanishing values of the 
induced dipoles p and m in the absence of external fields. 
This can only be accomplished if the denominator A± is 
zero in those equations, leading to 



(l/a E -G zz )(l/a M -G yy ) = H 2 



(22) 
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FIG. 9 (Color in online edition) (a) Lattice surface modes 
in a perforated semi-infinite perfect-conductor. The contour 
plot shows the modulus of the specular reflection coefficient 
for incident p-polarized light as a function of wavelength A and 
parallel momentum ku (see text insets for parameters). The 
upper-right inset shows a detail of the reflectivity as compared 
to the mode position predicted by Eqs. (23 1 and (24 1 (see 
arrow). A reflection coefficient larger than 1 is only possible 
for evanescent waves outside the light cone, (b) Lattice modes 
in a perforated thin film, as measured by Ulrich and Tacke 
(19721 (symbols). 



where we have set a' v = for infinitely deep holes (see 
Fig.[7|. The interaction sums G yy , G zz , and H are gener- 
ally small for s <C a, except near the lattice singularities 
discussed in Sec. |III.B[ In particular, near the light line 
for fcn > k, one has 



${G ZZ } « K{G yy } « ®{H} 



2nk 2 
k,a 2 



which corresponds to Eq. (14 1 with n — I — 0. Further 



more, upon inspection of an expansion for H similar to 
ufn, we find SH-ff} = outside the light cone, k« > k, and 



the remaining imaginary parts of all quantities in Eq. ( 22 1 
cancel out exactly because ^s{Gjj} = ^{a" 1 } = — 2/c T 73 
in that region. Combining these results, we obtain an ap- 
proximate long-wavelength dispersion relation from Eq. 



S 3 k 4 



ki = k 2 



r 



with 



4^2 
53" 



3?{l/a E } 3?{l/a M } 



(23) 



(24) 



Eq. (24) is exact in the s <§; a <C A limit, and it predicts 
the existence of lattice surface-bound modes under the 
condition 1/5R{1/q:e} + l/3?{l/aM} > 0. Here, we have 
used the area of the holes S to make T dimensionless. 



Calculated values of T arc offered in Fig. 10 'c) for var- 
ious hole geometries. The polarizability o?e («m) is ob- 
tained from the electrostatic (magnetostatic) far-field in- 
duced by an external electric (magnetic) field, as shown in 
Fig. 10 'a) [Fig. [lO^b)]. Interestingly, circular and square 
openings of the same area give rise to similar values of T. 
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FIG. 10 (Color in online edition) (a) Electrostatic electric- 
field flow lines for a circular hole drilled in a semi-infinite 
perfect-conductor subject to an external field E cxt perpendic- 
ular to the surface, giving rise to an electric dipole p = Ob S ext 
as seen from afar, (b) Magnetostatic magnetic-field flow lines 
for the same hole subject to an external parallel field H oxt 
and leading to a magnetic dipole m = qm H cxt . (c) Summary 
of polarizabilities for square and circular holes in perfect- 
conductor surfaces, normalized using the aperture area S. 
The values for the circular hole are taken from the h 2> b 
limit of Fig. [7| The circular opening in a thin screen is ana- 



lytical (Bethe 1944 Jackson 1999 1, but we must correct the 



right-hand side of Eq. ( 24 \ by a factor of 4 in this case because 



of cooperative interaction between both sides of the film. 



This parameter increases by an order of magnitude when 
the holes are made on thin screens instead of semi-infinite 
metals, producing lattice surface modes that are further 
apart from the light line (see Ulrich and Tacke 1972| , 
and therefore, more confined to the metal, as a result 
of cooperative interaction between both sides of the film 



[see a nalytical solutions for circular apertures ( Jackson 
1999) in last column of Fig. 10 'c)]. Another suggestive 



possibility is offered by split annular holes, which present 



resonant electric polarizability (Falcone et al 



2004), and 



by holes filled with high-permittivity materials (see Sec. 
III.E), for which the interaction with single-hole modes 
produces large departures of the extended surface states 
from the grazing light condition. 

Fig. [9^ a) shows calculated results for the reflection co- 
efficient of a drilled metal, obtained by rigorous solution 
of Maxwell's equations in which we use a plane-wave ex- 
pansion of the field outside the metal and a guided-mode 
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FIG. 11 (Color in online edition) Enhanced transmission 
driven by a localized resonance. The plot shows the normal- 
incidence transmission of a circular aperture drilled in a 
perfect-conductor film and filled with dielectric material for 
different values of the permittivity e (see labels). The trans- 
mitted power is normalized to the incoming flux within the 
hole area. 



modes are more bound in perforated thin films than in 
semi-infinite metals, as can be seen from the values of 
r given in Fig. 10 'c). Actually, the measured dispersion 



relation departs substantially from the light line close to 
the boundary of the first Brillouin zone. 



E. Interplay between lattice and site resonances 

The description of extraordinary optical transmission 
in terms of quasi-bound surface states driven by lattice 
singularities can be extended to other types of binding. 
In particular, a single hole filled with a dielectric of high 
permittivity can trap light in its interior, giving rise to 
cavity modes even for very subwavelength apertures, pro- 
vided the permittivity is sufficiently large to shrink the 
wavelength inside the dielectric to a value comparable 
to the diameter of the hole. This concept is explored 
in Fig. |11[ in which higher permittivities are seen to 
produce larger contraction of the wavelength inside the 
hole, so that the cavity mode condition is met at longer 
free-space wavelengths for fixed aperture size ( Garcia de 



expansion inside the holes (Garcia de Abajo and Saenz Abajo 2002 Garcia- Vidal et al. 2005). This process is 



2005). The lattice surface mode can be observed as a 



bright region with a dashed line showing the position at 
which the reflection coefficient becomes infinite. A de- 
tail of |r| for a specific wavelength (see dotted straight 
line) is shown in the inset. The position of the reso- 
nance predicted by Eqs. (23) and (p4| (see arrow in the 



accompanied by weaker coupling to external light (due 
in part to higher reflectivity of the dielectric-air inter- 
face), and therefore, narrower transmission resonances of 
increasingly larger height. Original predictions of this ef- 
fect (Garcia de Abajo 2002 1 have been recently corrobo- 



rated by experiment using microwaves (Garcia de Abajo 



inset) is in reasonably close agreement with the exact cal- 
culation, considering that the analytical model neglects 
neighboring-holes multipolar interaction, which is impor- 
tant for openings occupying 64% of the surface. Finally, 
Fig. [9jb) shows experimental results for a drilled thin 
film obtained by Ulrich and Tacke ( 1972 1. These surface 



et al. 2006 1 



An interesting situation is presented when localized 
modes like the ones just described are mixed with ex- 
tended lattice modes, like the surfac e states underly- 
ing extraordinary optical transmission (Garcia de Abajo 



et al. 


2006 


Ruan and Qiu 


2006 
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off in wavelength ( Jackson 1999 1. This led Roberts and 



McPhedran (1988) to theoretically explore the perfor- 
mance of periodic annular-hole arrays as band filters. 



j(a) p polarization 




More recently, Fan et al. ( 2005 1 have measured the in- 



k a fix 



creased transmission of infrared light assisted by these 
modes. Similar coupling to localized TEM modes occurs 
as well in slits, as we shall see in Sec. |III.F| 

The type of interplay phenomenon that we are describ- 
ing has been observed as well for localized and extended 
surface plasmons in the visible regime through the ab- 
sorption features of porous metals, in which Mie modes 
of spherical cavities in otherwise planar surfaces display 
a rich structure of hybridization and avoided crossings 



FIG. 12 (Color in online edition) Interplay between localized 
(site) and extended (lattice) resonances. The contour plots 
show the zero-order beam transmittance of a square array 
of circular holes drilled in a perfect-conductor film and filled 
with dielectric material of permittivity e = 50 as a function 
of parallel momentum kii and wavelength A. The orientation 
of k|| and the ratios between the hole radius b, the lattice 
constant a, and the film thickness h are specified in the in- 
sets. The light is p polarized in (a) and s polarized in (b). 
A transmission coefficient larger than 1 is only possible for 
evanescent waves below the light cone. 



both types of modes is illustrated in Fig. [12] through 
the zero-order transmittance of hole arrays filled with 
high-permittivity dielectric, calculated from the formal- 
ism presented in Sec. |III.C.3| All incident-light polar- 
izations interact with the cavity modes, giving rise to 
omnidirectional extraordinary transmission and invisibil- 
ity behavior near the individual hole resonance ( |Borisov 



(Baumberg 2006 Keif et al. 2005 2006 Teperik et al. 



2006ajb). The absorption can be even complete under at 



tainable experimental conditions (Teperik et al. 2005), 



implying black-body-like emission according to Kirch- 
hoff's laws of thermal radiation (Reif 19651. 



F. Slit and cylinder arrays 

Although we have extracted conclusions for particles 
and holes from his works, Wood reported his anomalies 
for ruled gratings rather than 2D structures ( |Wood 1902 
1935 1. 6 In fact, like gratings, cylinder and slit arrays ex- 



hibit lattice-resonance phenomena. But in contrast to 
holes, a single arbitrarily-narrow slit in a perfect conduc- 
tor supports at least one guided wave, the TEM mode 
( Jackson 1999 1 , which can couple to external p-polarized 



light (magnetic field parallel to the slit) giving rise to re 
cently predicted (Takakura 20011 and observed (Yang 



et al. 2005 Garcia de Abajo et al. 2005b). However, 
only p-polarized light couples to the n = 1, I = lat- 
tice singularity of Fig. g] which results in an avoided p/by lattice resonances fo/p pofarization (IPorto et al 



and Sambles 2002 ) Fabry-Perot resonances in transmis- 
As a consequence, light passage through slit ar- 
rays can be assisted either by coupling to the TEM mode 



avoided crossings have been recently found in microwave 
experiments ( |Hibbins et «?.[|2006[ ), confirming lattice sur- 
face modes and localized modes as two distinct mecha- 
nisms leading to enhanced transmission. 4 ' 5 Notice that 
s-polarized light is immune to the / — lattice singu- 
larities of Fig. [4] and this results in a reduced number 
of transmission features as compared to p polarization, 
in qualitative agreement with experimental observations 
(Barnes et al 



crossing of the hybridized modes [Fig. [12(a)]. Similar rigggl leading to similar interplay between localized anc 



extended resonances as discussed above (Marquicr et al. 



2005 ) . Incidentally, the analogy with annular hole arrays 



2004) 



Site resonances can occur in coaxial waveguides as well, 
via the so-called TEM mode, which does not have a cut- 



4 In a related context, avoided crossing of lattice modes are well- 
known to occur in coinciding Wood anomalies ( Stewart and Gall- 



is clear (see Sec. III.El. 

We shall consider first a periodic array of parallel nar- 
row cylinders, the axes of which define a single plane. 
Continuing with our tutorial approach, and focusing for 
simplicity on light incident with its electric field paral- 
lel to the cylinders, we note that Eqs. are still 
applicable here, provided qe and Q° are conveniently 
redefined. In particular, the polarizability has now di- 
mensions of area rather than volume, and it is given for 
instance by a E s = 7r6 2 (e— 1) for homogeneous cylinders of 
radius b and permittivity e ( |Bohren and Huffman 1983 1, 
with the optical theorem now leading to 3{l/aE| = 
—k 2 /A. The relevant dipole-dipole interaction compo- 



away|[T962} . 

incidentally, lattice modes are observed outside the light cone 
for p polarization. The transmission outside that cone is de- 
fined as the squared-amplitude ratio of incident and transmitted 
evanescent waves at the exit and entrance surfaces of the film, 
respectively. 



The reader is referred to the papers collected by |Maystre] ( |1993^ 
for an exciting historical overview of XX century milestones on 
gratings. 
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nent is given by the Green function of Helmholtz equa- 
tion in two dimensions, Q° = (ik 2 /tyH^^kR), where 
R is the distance measured in a plane per pendicular to 
the cylinders and H^ 1 is a Hankel function ( Abramowitz 
and Stegun, 1972). Then, proceeding with the lattice 
sum G(k||) in a way analogous to Eq. ([7b, one finds a 



relation similar to Eq. (Ill for the reflection coefficient 



of an array of lossless cylinders: 

-1 



1 + f*K{l/a E 



G(0)} 



Under normal incidence (fen =0), G is found to diverge 
as 

7T 1 

G(0) « , 

a 2 V2y/\/a-l 

for A > a, where a is the lattice period. This is similar to 
particle arrays [see Eq. ( 13 1] , so that the main conclusions 
from our previous discussion of those arrays apply here 
as well, and more precisely, the reflectivity can be made 
100% for arbitrarily narrow or weakly-scattering (e > 1) 
cylinders. 

A complete analysis along these lines has been recently 
reported for all possible incident polarizations ( Gomez 



IV. REAL METALS VERSUS PERFECT CONDUCTORS 

Metals of finite conductivity show significant differ- 
ences with respect to the perfect conductors consid- 
ered so far, the most remarkable of which is the exis- 
tence of intrinsic surface-plasmon excitations. The ba- 
sic understanding of these differences were laid out by 
Maystre| (|1972[) in the context of diffraction gratings (see 
also |McPhedran and Maystre||l974| and |Maystre[ 1984 1. 
Next, we shall examine (in a tutorial fashion) the con- 
sequences for the interaction between particles and holes 
decorating metal surfaces. 



A. Surface plasmons 

Conduction electrons in metals behave like a plasma 
that is capable of sustaining collective oscillations known 
as plasmons (e.g., longitudinal bulk modes, signalled 
by the vanishing of the dielectric function). The exis- 
tence of genuine surface plasmon oscillations was pre- 
dicted by |Rit chiej 1 1 95 7 [ and so on after confirmed by elec - 
tron energy- loss experiments (Powell and Swan 19591. 



Medina et al. 2006 Laroche et al. 



Since then, surface plasmons have developed into the 



rapidly growing field of plasmonics (Barnes et al. 2003 



similar lattice resonances, somewhat less pronounced, are 
obtained for E oxt perpendicular to the cylinders and with 
non- vanishing projection normal to the plane of the ar- 
ray. However, polarization components parallel to that 
plane and perpendicular to the cylinders cannot gener- 
ate lattice resonances, because the interaction between 
distant dipoles aligned with their separation vector R. 
decays as 1/i? 3 / 2 in 2D, which is insufficient to produce 
a divergence in G. 7 

Finally, we can establish a relation between cylinder 
arrays and slit arrays using arguments similar to those of 
Sec. |III.C. 1| for particle and hole arrays. More precisely, 
a slit array cut into a thin metal screen and illuminated 
with E cxt perpendicular to the apertures can be analyzed 
using the above results as applied to the Babinet-related 
stripe array (i.e., a periodic array of stripes laying on a 
single imaginary plane) for E cxt parallel to the stripes. 
Under normal incidence, the required component of the 
polarizability reads a E w — 27r/fc 2 [ln(fcfe/8) + 7 + in/2], 
where 7 = 0.57721 is the Euler constant and b <C A is 
the stripe width (van de Hulst[ 1981). Interestingly, a E 



2006J), suggesting that pzbayl [20061 |Zia et al] [2006[) owing to their p otentia 



applicabili t y to areas as diverse as biosensing ( |Schus-| 
ter et al] 1993|), signal processin g through plasmonic 



circuits (|Bozhevolnyi et ~aT. 2006), or laser technology 



(Colombelli et al. 20031 



Planar surfaces possess translational invariance that 
provide plasmons with well-defined parallel momentum 
fen exceeding that of light outside the metal and thus 
becoming truly surface-bound modes. Their dispersion 
relation can be readily derived from the divergence of 
the Fresnel coefficients for p polarization (surface-bound 
fields without external sources), leading to (Raether 



1988) 



uSP 



1 



(25) 



for a metal-air interface. This surface plasmon dispersion 



relation is represented in Fig. 13 'a) for a Drude metal de 



scribed by Eq. (12). In the long fcy limit, the surface plas- 
mon frequency satu rates to Ritchie's non-retarded plas- 
mon (iRitchiel 119571). 



diverges in the electrostatic limit, so that even a single 
narrowing slit will exhibit a divergent transmission cross 
section. This scenario can be traced back to the above- 
mentioned site resonances produced by the TEM mode of 
slits in thick screens. As a consequence, the interaction 
between slits can be very large, resulting in strong red 
shifts of the transmission peaks relative to the Rayleigh 
condition. 



Surface plasmons are characterized by three different 
length scales, as depicted in Fig. [BJb): their propa- 
gation distance along the surface (~ l/2!5{fc| p }), their 
penetration into the surrounding medium (~ 1/23{A:_l}, 
where k± = — fc/ye + 1 is the normal momentum), 
and their penetration into the metal (the skin depth 
~ l/23{— efcj^}). Interestingly, the interaction between 
plasmons in either sides of a thin film gives rise to two 
plasmon branches, as measured by electron microscopy 
(Pettit et al. 1975 Vincent and Silcox 19731, one of 



7 This is because X^°° 1/n 3 / 2 is finite. See also Sec. III.B 



which has been found to propagate along very long dis- 
tances thanks to exclusion of the electric field from the 
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FIG. 13 (Color in online edition) (a) Surface plasmon dispersion relation for a Drude metal of bulk plasmon frequency ui p . 
(b) Extension of the plasmon field into the metal (skin depth), into the vacuum, and along the surface (propagation distance) 
for several metals, as obtained from measured optical constants (Johnson and Christy 1972 Palik 19851. 



metal (JSarid 
have 



1981[ ). Well defined plasmons require to 
<C 3?{— e}, but similar long-range surface- 



exciton polaritons exist in thin films for 3{e} |5i{e}| 



(Yang et al. 19901 



Features in metal surfaces produce scattering of plas- 
mons in a similar way as light is dispersed by particles. 
This is actually a way to couple externally incident light 



to plasmons, for instance using gratings ( jLoewen et al. 



1984 Ritchie et al. 19681. We find a neat demonstra- 



tion of these ideas in the observation of surface-plasmon 



bands for periodic surface decoration (Kitson et al. 


1996 


Ritchie et al. 


1968 Stewart and Gallaway 


1962 


1 and 



in the reflection of surface plas mons at point scatter - 
ers arranged as parabolic mirrors (iNomura et al. 20051. 



Similarly, holes perforating films have strong influence 
on surface plasmons, which play an important role in 

How- 



their optical transmission (Ghaemi et al. 
ever , in the perfect-conductor limit 
k 



with 



1998) 



Eq. 



(251 yields fc| p = k, with zero skin depth and infinite 
penetration into the vacuum, that is, there are no longer 
surface-bound modes. In the following we shall explore 
the transition between plasmonic and perfect-conductor 
regimes, in an attempt to clarify seemingly contradictory 
statements regarding the role of surface plasmons to en 
hance flSchroter and Heitmann 1998 1 or to suppress ( Cao 



and Lalanne 2002 ) extraordinary optical transmission in 
striped thin films, or the heated debate opened by the 
explanation of recent outstanding experiments dealing 
with the interaction between a slit and a groove ( Garcia 



Vidal et q?.||2"006]|Gay eFaT|[2006l[Lalanne and Hu gonin 
20061). 



stance can lead to a net parallel electric dipolc in a thin 
metallic film (Rindzcvicius et al. 20071. 



In order to illustrate this concept, we have considered 
in Fig. [14] the effective polarizability of a silver spherical 
particle in front of a silver surface for a constant ratio of 
the radius to the wavelength, b/X = 0.1. We can observe 
an electric Mie mode (Mie 19081 in the visible, accom- 
panied by negligible magnetic response. However, the 
metal behaves increasingly closer to a perfect conductor 
at longer wavelengths, so that currents compete eventu- 
ally with polarization, thus displaying magnetic polariz- 
ability that becomes cvm = —b 3 /2 for an isolated perfect- 



conductor sphere in the long- wavelength limit (Jackson 



1999), to be compared with the electric polarizability 
«e = b 3 . Nevertheless, the latter is quenched by prox- 
imity of the metal flat surface under normal-incidence 
illumination conditions. The onset of magnetic response 
occurs when the particle becomes large compared to the 
skin depth ~ 20 nm [see Fig. 13 'b)]. These results fol- 
low from dipolar Mie scattering, conveniently corrected 
by surface reflection coefficients, which qualitatively de- 
scribe the polarizability strength of the coupled particle- 
surface system. 

This has important consequences for understanding 
patterned surfaces and hole arrays. Electric dipoles dom- 
inate the response of features smaller than the skin depth, 
whereas magnetic dipoles can be significant for larger 
sizes, and only parallel electric dipoles and perpendic- 
ular magnetic dipoles survive in the limit of negligible 
skin depth. We are of course restricting our discussion 
to particles or apertures that are small compared to the 
wavelength, but these conclusions can be generalized to 
higher-order multipolcs for bigger features. 



B. Polarization schemes 



The condition that parallel electric dipoles and per- 
pendicular magnetic dipoles are excluded from perfect- 
conductor surfaces (see Fig. 10 1 is relaxed in metals of 



finite conductivity. Polarization charges in a hole for in- 



C. Dipole-dipole interaction 

New dipole orientations and the presence of surface 
plasmons in real metals demand that we revisit the in- 
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FIG. 14 (Color in online edition) Effective polarization 
strength of a silver sphere near a silver planar surface. The 
sphere radius is a tenth of the wavelength. The polarization 
is normalized to the sphere volume. The dielectric function 
of silver is taken from [Joh nson and Christy, |1972| 



teraction between features in tailored surfaces. In par- 
ticular, the dipolar field in free space, which decays away 
from the source as 



g° 



rdkR 



R 



(26) 



and governs the interaction between small features in 



perfect-conductor surfaces (see Sec. III. A I, must be sup- 
plemented by reflected fields near real metals, leading to 
an interaction tensor of the form 

g = g° + g r . 

As a result, light impinging on a hole can couple to cir- 



cular surfacc-plasmon waves (Chang et al.\ 


2005 Popov 


et al. 


2005 


Wannemacher 


2001 Yin et al.\ 


2004), whose 



field strength shows a rather different decay dependence 
with distance as 



JfcPfl 



Vr 



(27) 



This expression is consistent with energy flux conserva- 
tion for any surface-bound mode, 8 with dissipation de- 
scribed through the imaginary part of fc? p . The slow 



drop of Eq. (27) with distance compared to Eq. (26 1 can 



explain the observed enhancement of the interaction be- 
tween small particles in plasmonic metals ( |Stuart and 
Hall[ |1998 |, and it is illustrated in Fig. 15 showing the 



field produced by a dipole near a metallic surface as cal- 
culated from a trivial extension of our tutorial approach 
formalism presented below. 

The interaction between pairs of electric and magnetic 
dipoles near a metal surface is analyzed in detail in Fig. 
16 'a) for all possible orientations except perpendicular 
magnetic dipoles, which are forbidden in perfect conduc- 
tors and should take small values in real metals. More- 
over, symmetry forbids the interaction of all other pairs 
that are not shown in the figure. For surface features in- 
ducing electric dipoles under normal incidence in a plas- 
monic metal (see Fig. 14), the dominant interactions orig- 
inate in electric-dipole pairs aligned with their separation 
vector R (see Fig. 



16) 



quite different from perfect con- 
ductors, which are governed by magnetic dipoles perpen- 
dicular to R. However, the latter can contribute in plas- 
monic materials as well for large features compared to 



the skin depth, as we discussed in Sec. IV. B As a thumb 
rule, the mutual dipole orientations that lead to the long- 
range interaction dependence given by ( 27 ) are compati- 



ble with non- vanishing surface-plasmon field components 
emanating from those dipoles [i.e., plasmons with m = 
azimuthal symmetry for normal electric dipoles, like in 
Fig. 15 or m = ±1 for parallel dipoles]. 



The interaction between dipoles in front of a planar 
surface admits a representation in parallel momentum 
space similar to Eq. |5|, but involving now the Fresnel 
reflection coefficients for s and p p olarization ( Blanco and 
Garcia de Abajo] [2004] [Weyi] p!9| , r s = (fc z - fc'J/(fc z T 
k' z ) a nd r v = (efc z — fc'J /(ek z + k' z ), respectively ( Jackson 
1999| ), where k z = y/k 2 - Q 2 and k' z = y/k 2 e-Q 2 . In 
particular, for elect ric dipoles parallel to the surface x 
direction, one finds ( |Ford and Weber 



1984 Weyl 1919) 



g r xx 



d 2 Q 

2nJ k z Q 2 



i(Q-R+fe i | 2: |)j fc 2g2 



k z Q x 



(28) 



where z is the sum of distances from the dipoles to the 
surface, and we are interested in the z — > limit. This 
expression is general and leads to g xx = in perfect 
conductors, for which r p = — r s = 1. 

The strong surface-plasmon-mediated interaction de- 
scribed by Eq. (27 1 arises from the pole of the Fresnel 
= fcn P , which admits the Laurent ex- 



coefficient r p at 



pansion (Ford and Weber 1984) 



2Bk 

Q-k 



SP ' 



(29) 



8 The Poynting vector produced by a dipole when the fields are 
propagated by means of Eq. (127b dies off as 1/R, if we neglect 
the attenuation produced by Sjfc^ p }. Then, the integral of the 
radial Poynting vector over a circle of radius R centered around 
the dipole and lying on the surface is independent of R, so that 
the photon flux is conserved, indicating that we are dealing with 
surface-bound propagation. 



with 



B = [e/(l + e)f/ 2 /(l-e). 



Performing asymptotic analysis for large R and retaining 
only the contribution from this pole in the integral of Eq. 
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FIG. 15 (Color in online edition) Instantaneous electric field set up by a perpendicular electric dipole (see vertical arrows) 
sitting at distance A/20 from the surface of a metal described by Eq. ( 12 1 with uj p = 15 eV and damping r\ = 0.6 eV (typical of 
Al) at frequency u = cu p /2. The electric-field component parallel to the surface (this is radial with respect to the position of the 
dipole) and the component along the surface normal are represented separately. Poynting vector flow lines are superimposed 
on the plot of the normal component. 
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FIG. 16 (Color in online edition) (a) Schematic represen- 
tation of the scaling of dipole-dipole interactions for electric 
and magnetic dipoles with respect to their separation R near a 
metallic surface. The interaction decays as exp(ifci?)/_R™ near 
a perfect conductor or as exp(ifcy P R) / R m near a metal with a 
dominant surface plasmon (see text insets for values of the ex- 
ponents n and m). (b) Dipole-dipole interaction near a silver 
surface at a wavelength of 750 nm (three upper solid curves) 
as compared with the plasmon-pole approximation (three up- 
per dashed curves, see text). We also show the interaction at 
a wavelength of 10 mm (lower curve, perfect-conductor limit). 
The dipole-dipole separation vector R is taken along x. 



p8| (plasmon-pole approximation; see Ford and Weber 
19&*l, we obtain 9 



e+1 



H^\kfR) + H^\kfR) { ^ 1 ^l 



2Txk A B 2e e 



e + 1 V \-nk 



SP 



R 



(30) 



where the second approximation comes from the asymp- 
totic behavior of Hankel functions for large arguments 
(Abramowitz and Stegun 19721, so that one obtains the 
result anticipated in Eq. (|27p. The above approximate 
expression in terms of Hankel functions is compared with 
the direct numerical evaluation of Eq. ( 28 1 , and similar 

m Fig. [lib). 



expressions for other dipole orientations, 
The agreement at A = 750 nm is excellent for R > A, in- 
dicating that lattice resonances in an array will be really 
dominated by surface plasmons at that wavelength. Fig. 
16 illustrates as well a much faster decay of Q yy as 1/ i? 3 / 2 



Tot electric dipoles oriented orthogonal to R and parallel 
to the surface, and as 1/R for normal electric dipoles in 
the perfect-conductor limit. 



D. Discrepancies in lattice resonances and enhanced 
transmission 

The dissimilar behavior of plasmonic metals and per- 
fect conductors discussed in the previous sections leads 
to qualitative differences in extraordinary optical trans- 
mission, arising in part from the 1/(Q — fc| P ) dominant 
pole of the inter-hole interaction in momentum space [see 



Eqs. (28l-(30 1] 



Considering for simplicity a square array under normal 
incidence, we can analyze the lattice sum in a real metal 



9 It should be noted that the asymptotic behavior of Q° [see Eq. 
( |26| ] comes from the Q = k region of the integral in Eq. (JsJ 
and responds to the pole l/fc z . This pole is canceled exactly by 
Eq. n28ll , in which r p = r s = — 1 at grazing incidence (i.e., for 
Q = k). Therefore, the only relevant contribution to Q for large 
R originates in the plasmon pole of Q r . 
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[i.e., Eq. Q with Q substituted for Q°] following the pro- 
cedure that led to Eq. Q, but starting now from Eqs. 
^ and (28 1. In a diffrationless array, there are just two 
identical singular terms in the corresponding sum over 
reciprocal lattice vectors, leading to 



G 



EE 

xx 



c 



SP 



47ry A 
aX J Asp /a 



1 



(31) 



> 



for Jft{A SP } 

plasmon wavelength and C 



where Agp = 27r//c^ p is the surface- 

We have ex- 



LB/Ve+1 

plicitly indicated with superscripts that describes 
the interaction between electric dipoles (E), which can 
coexist with parallel magnetic dipoles (M) (see Fig. 16 1. 
The remaining relevant lattice sums are G 

1) and G® y = —Gys ^ ~ V e + 1 G 



EE 



iMM 

T yy ~ v c "t 
Now, the for- 



malism presented in Sec. III. A can be easily extended to 
patterned surfaces and hole arrays in real metals using 
these expressions of the lattice sums rather than those 
for perfect conductors. 10 

In the polaritonic regime of surface plasmons, in which 
their dispersion relation approaches the light line (see 
e| is large and the dominant lattice sum scales 
j 1 in the plasmon-pole approximation, so 

that for sufficiently high |e| the perfect-conductor limit of 



Fig 
as G 



Eq. (13) dominates over the plasmon. 



A descriptive example of the transition from plasmonic 
to perfect-conductor behavior is offered in Fig. 
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which the energy released by a dipole sitting near a sur- 
face is divided into plasmon launching (Isp) and emit- 
ted light (/free)- This relates to the question, which of 
the two mechanisms (plasmons or propagating radiation) 
produces stronger interaction with a nearby surface fea- 
ture. Plasmon launching dominates near the electrostatic 
plasmon, reaching an efficiency close to 100% in silver. 
As the wavelength advances towards to infrared, the plas- 
mon is less bound to the surface and has weaker coupling 
to our dipole. As an example of application, when light 
pops out of a narrow hole after being guided through a 
TE mode (e.g., in a circular hole infiltrated with a di- 
electric of refraction index n ^> 1 and for X/n < 3.46), 
the equivalent dipole describing the hole lies parallel to 
the surface. That is the situation depicted in the inset of 
Fig. [17} 

A more explicit comparison of discrepancies between 
both metallic regimes for holes is offered in Fig. [18] which 
shows the lattice sum for parallel magnetic dipoles [ob- 
tained by summing Eq. (28) for gold, with the expres- 
sion in square brackets replaced by k 2 Q 2 x r p — k 2 Q 2 r s ], 
together with a geometrical construction like in Fig. [6] 
applied now to two different aperture sizes. It should 



10 Our analysis can be applied to metals embedded in a dielectric 
of refraction index n simply by using the reduced wavelength 
X/n everywhere instead of A and by interpreting e as the ratio of 
permittivities in the metal and in the dielectric. 
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FIG. 17 (Color in online edition) Relation between the power 
radiated after transmission through a deep subwavelength 
hole (/free) and the power emanating as surface plasmons 
(Isp) for gold and silver, derived in the small-hole limit. The 
metal dielectric function has been taken from IJohnson andl 
Christy! il972 I. 



be noted that the exact calculation (solid curves) com- 
pares extremely well with analytical expressions [sym- 
bols, obtained from Eq. ( 13 1 for the perfect conductor 



and from Eq. (31 1 for the plasmonic metal, which needs 
to be multiplied by — (e+ 1) in order to apply it to mag- 
netic rather than electric dipoles]. The lattice sum sin- 
gularity in perforated gold takes place to the red as com- 
pared to the perfect-conductor case, because the surface- 
plasmon wavelength is shorter than the light wavelength 
in the surrounding dielectric. Moreover, the lattice sum 
diverges as X/^J X/n — a and 1/(Asp — a) in perfect con- 
ductors and plasmonic metals, respectively, according to 
Eqs. (13) and (pi]), thus leading to different dependence 



of the position of the lattice surface resonance on hole 
size (see points of intersection with horizontal lines in 
Fig. 18 1; the lattice resonance is further away from the 



interaction sum singularity (and a given change in hole 
diameter produces larger peak shift) in the plasmonic 
case considered in the figure. 

The crossover between both types of behavior is ex- 
plored in Fig. 19 through the absorbance of (i) a silver- 
particle array in silica, (ii) the same array near a silver- 
silica interface, and (iii) an array of silica inclusions right 
underneath the metal-dielectric interface. We have done 
these calculations using a layer KKR method to solve 
Maxwell's equations (|Stefanou et al.\ |1998l 120001). In 



the case (i) a maximum in absorption occurs near the 
Rayleigh condition for light propagating in silica (i.e., 
X/n = a), whereas case (iii) shows a single maximum 
shifted to the right of the Rayleigh condition for the pla- 
nar interface plasmon (Asp = a) (Ghaemi et al. 19981. 
The conclusion is that plasmons are mediating the inter- 
action among the dielectric inclusions, with no signature 
of any anomaly near X/n = a whatsoever. An intermedi- 
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FIG. 18 (Color in online edition) Lattice sums and lattice 
resonances in a square array of holes drilled in gold vs a per- 
fect conductor. The real part of the exact lattice sum for 
interaction of parallel magnetic (M) dipoles is shown for gold 
(black curve) and for a perfect conductor (PC, grey curve), 
as compared to analytical approximate expressions (symbols). 
The Rayleigh condition for a period a = 600 nm is indicated 
by black and grey vertical dashed lines for light in the dielec- 
tric (A/n = a) and for surface plasmons (Asp = a), respec- 
tively. Changes in the inverse magnetic polarizability of cir- 
cular holes of different size [horizontal lines, as obtained from 
Fig. [TJb)] , lead to different wavelengths of the lattice surface 
modes, as indicated by vertical arrows for the condition that 
the real part of the denominator of Eq. (|3| be zero. 



ate situation is encountered in case (ii), showing features 
near the two types of Rayleigh conditions. 

It should be noted that Asp has an imaginary part aris- 
ing from absorption, and although it is small for noble 
metals, in which plasmons can travel long distances along 
the surface, as shown in Fig. 13 ^b), we find that Eq. (31 ) 
does not describe a divergence, but rather a Lorentzian 
of finite width. This affects the height of the transmis- 
sion maxima, below 100% in lossy metals. Furthermore, 
apertures perforated in metals of finite conductivity will 
appear to be wider by the skin depth effect, and their 
effective polarizability must be lossy. 

Without entering into further considerations regarding 
how finite conductivity affects the hole polarizability, let 
us just point out that the wavelength at which the noted 
intersection takes place in Fig. 18 (i.e., the wavelength 
of the lattice surface-bound mode) is in excellent agree- 
ment with the transmission peaks measured by Krishnan 
et al. (20011 and reproduced in Fig. [l|a). The verti- 
cal arrows in that figure indicate the predicted positions 
of the transmission maxima, obtained by increasing the 
hole size by the skin depth to an effective diameter of 
250 nm. This agreement is remarkable, given our neglect 
of higher-order multipolar terms in the hole polarization. 
The shift with respect to the Rayleigh condition for sur- 
face plasmons (vertical solid lines in Fig. [lj is signifi- 
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FIG. 19 (Color in online edition) Normal-incidence ab- 
sorbance of (i) a silver particle array embedded in silica (re- 
fraction index n = 1.45), (ii) the same array near a planar 
silver-silica interface, and (iii) an array of silica inclusions 
buried in silver below a silver-silica interface. All particles 
are spheres of 200 nm in diameter. The arrays have square 
symmetry with lattice constant a — 500 nm. The distance 
from the sphere surfaces to the planar interface is 10 nm in 
the buried silica particles and 900 nm for the silver particles. 
The Rayleigh conditions for the reduced wavelength of light 
in the silica (A/n = a) and for the wavelength of the silver- 
silica interface plasmon (Asp = a) are indicated by arrows A 
and B, respectively. 



cant, triggered by large, plasmon- mediated interaction 
between apertures, as explained above. Similar conclu- 
sions can be drawn for the silver film of Fig.[ljb), in which 
the results from the above analytical model are shown as 
dashed curves (divided by a factor of 5). Only magnetic 
dipoles are taken into account, with the hole polarizabil- 
ity calculated for a perfect conductor. The transmittance 
is obtained from Eq. (l9| with G xx replaced by its plas- 
monic counterpart, G y l y ■ Although the Rayleigh con- 
dition for plasmons (solid vertical lines in Fig. [T| agrees 
only with the transmission minima in silver (presumably 
because gold is more dissipative in this spectral region, 
so that the polarizability of the holes requires a more re- 
alistic description including absorption), the comparison 
with experiment is excellent, given the simplicity of the 
analytical model, which should become exact in the limit 
of small scattering features (e.g., for nanoparticle arrays 
on a metal substrate). 



V. CONCLUSION 

Light scattering in planar periodic systems gives rise to 
resonant phenomena that have common origins in parti- 
cle and hole arrays, both for reflection and for transmis- 
sion. Namely, (i) the interaction between lattice sites 
shows a divergent behavior when a diffracted beam be- 
comes grazing ( |Lord Rayleigh 1907), producing a min- 
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imum in both the reflectivity of particle arrays and the 
transmission of hole arrays; (ii) a lattice resonance can 
be established at a wavelength to the red of that condi- 
tion, leading to maxima in both the reflectivity of par- 
ticle arrays and the transmission of hole arrays; (iii) 
these effects have the same origin as Wood's anoma- 



lies (Wood 1935) and th ey can be d escribed in the lan- 
guage of Fano lineshapes (Fano 1961 ); (iv) the noted lat- 



tice resonance persists for incident evanescent light, with 
the reflectivity's becoming infinite in non-dissipative sys- 
tems (e.g., patterned perfect-conductors, but also pat- 
terned dielectrics) , thus defining truly surface-bound 



states ( Garcia de Abajo and Saenz 


2005 


Hibbins et al. 


2005 


Pendry et al.\ 2004 


Ulrich and Tacke 


19721; (v) 



these extended lattice resonances mix strongly with other 
modes localized at specific sites, like those created by 



nanoparticle and nanovoid plasmons (Keif et al. 2006 
|Teperik et al. 2006a); (vi) for metals with well-defined 
surface plasmons, the interaction between holes or par- 
ticles in the vicinity of the surface is mediated by these 
excitations, so that we have to reformulate the condi- 
tion of a diffracted beam's becoming grazing using the 
surface plasmon wavelength rather than the incoming or 
transmitted light wavelength. 

We have shown that particle arrays and hole patterns 
in perfect conductors share in common the asymptotic 
form of their interaction, summarized by Eq. (26), which 



produces singularities at the Rayleigh condition when 
summed over the lattice, for instance for A — a under 
normal incidence on square arrays, and gives rise to sur- 
face states at slightly larger wavelengths. However, the 
plasmon-mcdiated interaction in noble metals is more in- 
tense, as shown in Eq. (27), thus producing sharper diver- 
gences and stronger collective interaction. In this case, 
the singularities occur at the band-folded plasmon lines 
(e.g., when Asp = a under normal incidence on square 
arrays), and the lattice surface-bound states (i.e., the 
plasmons of the patterned metal) exist again to the red 
with respect to those lines. 

All of these effects have been described here within 
a common tutorial approach based upon interacting 
dipoles that is not only able to explain the observed ef- 
fects; its simplicity has allowed us to extract some sur- 
prising conclusions. One of them is that arbitrarily- weak 
scatterers forming a periodic structure and made of non- 
dissipative materials can also produce intense lattice res- 
onances: given an array of arbitrarily-small particles of 
positive polarizability, it is always possible to find a wave- 
length (close to the period for square symmetry and nor- 
mal incidence) at which light is totally reflected; accord- 
ingly, it is possible to obtain full transmission through 
holes however narrow, drilled in arbitrarily-thick perfect- 
conductor films. 

Interestingly, the lattice periodicity alone determines 
the magnitude of the induced dipoles needed to pro- 
duce complete reflection by small particles or total trans- 
mission through narrow holes. Moreover, the polariz- 
ability scales with the cube of the hole/particle diame- 



ter. Combining these two statements, we find that the 
self-consistent electric field acting on particles or aper- 
tures under such resonant conditions increases when they 
shrink and can reach extremely high values only limited 
by absorption and lattice imperfections, thus opening 
new possibilities for applications in nonlinear all-optical 
switching and biosensing. 

The simplicity and power of the model that has been 
presented here will surely find application to explain 
many other effects related to light scattering in planar 
periodic systems and can be inspiring for devising new 
phenomena. 
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